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ABSTRACT 

The accumulation of Swift observed gamma-ray bursts (GRBs) gradually makes it 
possible to directly derive a GRB luminosity function (LF) from observational lu- 
minosity distribution, where however two complexities must be involved as (i) the 
evolving connection between GRB rate and cosmic star formation rate and (ii) ob- 
servational selection effects due to telescope thresholds and redshift measurements. 
With a phenomenological investigation on these two complexities, we constrain and 
discriminate two popular competitive LF models (i.e., broke-power- law LF and single- 
powcr-law LF with an exponential cutoff at low luminosities). As a result, we find 
that the broken-power-law LF could be more favored by the observation, with a break 
luminosity Lt — 2.5 x lO^^erg s^^ and prior- and post-break indices i^i ~ 1.72 and 
P2 = 1.98. For an extra evolution effect expressed by a factor (1 -I- 2)*, if the matallic- 
ity of GRB progenitors is lower than ^ O.IZq as expected by some coUapsar models, 
then there may be no extra evolution effect other than the metallicity evolution (i.e., 
S approaches to be zero). Alternatively, if we remove the theoretical metallicity re- 
quirement, then a relationship between the degenerate parameters S and .^max can be 
found, very roughly, 6 ~ 2A{Z^i^x/Zq — 0.06). This indicates that an extra evolution 
could become necessary for relatively high metallicities. 
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1 INTRODUCTION 

Some confirmed associations between gamma-ray burstj^ 
(GRBs) and Type Ib/c supernovae (Stanek et al. 2003; 
Hjorth et al. 2003; Chornock et al 2010) robustly suggest 
that GRBs are powered by the collapse of the core of mas- 
sive stars, which is also widely accepted in theory (Woosley 
1993; Paczyhski 1998; Fryer et al. 1999; Wheeler et al 2000; 
Woosley & Bloom 2006). In other words, the detection of 
each GRB provides a witness of the death of a massive star. 
Moreover, the intense brightness of GRBs makes them de- 
tectable even at the edge of the universe (the highest redshift 
of GRBs is ~ 9.4 as reported by Cucchiara et al. 2011). 
So GRBs can in principle be used as a tracer of the cos- 
mic star formation history. The crucial problem is whether 
GRBs are an unbiased tracer or, more directly, how to cali- 
brate GRB event rate to star formation rate (SFR). On one 
hand, the cosmic evolution of metallicity could be involved. 
This is because a very high angular momentum is required 
for GRB progenitors and meanwhile massive stars in lower- 
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^ Throughout we refer only "long" gamma-ray bursts with T90 > 
2 s, where Tgo is the interval observed to contain 90% of the 
prompt emission. 



metallicity environments are less likely to loss much angular 
momentum through stellar winds (e.g., Meynet et al. 1994; 
Langer & Henkel 1995; Vink & de Koter 2005; MacFadyen 
& Woosley 1999; Wooseley & Heger 2006). On the other 
hand, the luminosity function (LF) of GRBs can also play 
an important role in the conversion from the observed GRB 
redshift distribution to GRB formation history, since the lu- 
minosity selection by telescopes leads to a lower detection 
probability for higher-redshift GRBs. 

To directly derive a GRB LF was impossible before the 
launch of Swift (Gehrels et al. 2004), since there were only 
quite a few GRBs whose redshifts are measured. A possi- 
ble alternative way invokes some luminosity-indicator re- 
lationships to avoid redshift measurements (e.g. Yonetoku 
et al. 2004; Fenimore & Ramirez-Ruiz 2000; Firmani et al. 
2004), but the robustness of those indicators may be not 
high enough. A much more popular method is to assume a 
LF form with a few model parameters and then to fit the flux 
distribution of the observed GRBs (log A'' — log P distribu- 
tion; Schmidt 1999; Porciani & Madau 2001; Firmani et al. 
2004; Guetta et al. 2005; Natarajan et al. 2005; Daigne et al. 
2006; Salvaterra & Chincarini 2007; Salvaterra et al. 2009; 
Campisi et al. 2010). Correspondingly, before Swift, a very 
large sample of GRBs had been provided by the Burst and 
Transient Source Experiment (BATSE) on board Compton 
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Observatory. Nevertheless, by such a fitting to log A'^ — log P 
distribution only, it is not easy to eliminate the degener- 
acy among the model parameters and even to determine the 
form of the LF. As a result, two competitive LF models as 
a broken-power law (BPL) and a single-power law with an 
exponential cutoff (SPLEC) at low luminosities are usually 
adopted in literature. 

Thanks to Swift spacecraft, in the past few years the 
number of GRBs with measured redshift grows rapidly. This 
makes it possible to provide more stringent constraints on 
the LF parameters (Daigne et al. 2006; Salvaterra & Chin- 
carini 2007; Salvaterra et al. 2009; Campisi et al. 2010). The 
new constraints robustly rule out the models in which GRBs 
unbiased trace the cosmic star formation or GRBs are char- 
acterized by a constant LF. In other words, an evolution 
effect is suggested. Li view of the not-small size of the Swift 
GRB sample, it has become possible to derive a GRB LF 
only with Swift GRBs. Very recently, Wanderman & Piran 
(2010) tried to directly convert the luminosity distribution 
of Swift GRBs to a LF, without a prior assumed LF form 
and without a help from the BATSE data In such a LF- 
determination process, the treatments on observational se- 
lection effects play a very important role. Meanwhile, a pos- 
sible extra evolution effect should also be paid much atten- 
tion to. In this paper, with a phenomenological investigation 
on the evolution effect and the selection effects, we constrain 
and discriminate the BPL and SPLEC models by using Swift 
observed GRBs. 

In the next section, some observational and theoreti- 
cal materials related to the GRB LF are described. In sec- 
tion 3, the evolution effect is constrained and analyzed with 
relatively-high-luminosity GRBs. In Section 4, firstly, we 
derive an initial LF in both the BPL and SPLEC models 
by directly fitting the observational luminosity distribution 
of GRBs. Secondly, we analyze and constrain the so-called 
redshift-desert effect with the initial LFs. Finally, the com- 
bination of the above two processes gives a final constraint 
on the GRB LF, with which the prior selected luminosity 
threshold is checked. In Section 5, conclusion and discussion 
are given. 



2 BASIC MATERIALS 
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Figure 1. Luminosity-redshift distribution of 125 Swift GRBs 
with redshift, where the shaded region represents the luminosity 
threshold adopted in our calculations (see Equation |9] and ex- 
planations there). The 24 data labeled by open asterisks will be 
excluded from our most statistics. 
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2.1 Swift observed GRBs 

In the past few years. Swift has greatly promoted our 
understanding of GRBs. Here we take GRBs with mea- 
sured redshift z from the Swift archiv(j3- For most of these 
GRBs till GRB 090813, their spectral peak energy Ep and 
isotropically-equivalent energy release _Biso in the burst rest- 
frame 1 — 10* keV band have been provided by Butler et al. 
(2007, 2010). However, it should be noticed that, due to the 
narrow energy bandpass of the Swift Burst Alert Telescope 
(BAT), the burst spectral parameters in Butler et al. (2007, 
2010) are actually estimated by the Bayesian statistics, but 
not directly observed. Anyway, as in Kistler et al. (2008, 
2009) and Wang & Dai (2009), an average luminosity can 
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Figure 2. Upper panel: Comparison between the distributions 
of average fluxes of the total 406 GRBs and the 125 GRBs with 
redshift; Lower panel: Probability of redshift measurements as a 
function of flux and an empirical fitting given by Equation ((2]|, 
where the horizon error bars correspond to the bin width. The 
vertical dashed line represents the selected threshold = 5 X 
10^^ erg s^-'^cm^^. 



be roughly estimated for these GRBs by 



E: 



^[l-lO^lkcV] 



so[l-10lkGV] 



no /{I 



(1) 



In our statistics, GRBs with L < 10 erg s^ will be ex- 
cluded, since they may belong to a distinct population called 
low-luminosity GRBs (Soderberg et al. 2004; Cobb et al. 
2006; Chapman et al. 2007; Liang et al. 2007). Finally, 125 
GRBs are selected and their luminosity-redshift distribution 
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is shown in Figure 1. Correspondingly, including the GRBs 
without redshift, there are totally 406 GRBs detected by 
Swift till GRB 090813. 

In the upper panel of figure 2, we present the distribu- 
tions of the observed average 

i5-i50kcv], for both 

the total 406 GRBs and the 125 GRBs with redshift. The 
ratio between these two distributions generally displays the 
flux-dependence of the GRB redshift measurements and, as 
shown in the lower panel, such a redshift detection proba- 
bility can be empirically expressed by 



•dziP) = min 



0.27 + 



2.0 X 10-^erg s^^cm^^ 



(2) 



A similar result has also been given in Qin et al. (2010). 
On the other hand, the redshift detection probability may 
also depend on redshift itself, which will be investigated in 
section 4.2. 



2.2 Model 

The luminosity-redshift distribution of GRBs is determined 
by both the LF </'(L) and the comoving rate R{z) of GRBs, 
which are respectively deflned by 



(3) 



(4) 



dL ='^(^)' 
and 

dN _ ■ dV{z)/dz 
az 1 + 2 

where the dot represents time derivation, the factor (1 + z) 
is due to the cosmological time dilation of the observed 
rate and dV{z)/dz is the comoving volume element. In 
the standard A-cold dark matter cosmology, dV{z)/dz — 
ATvdc{z)'^c/ H{z) with dc{z) — di{z)/{l + z), where the lu- 
minosity distance reads di{z) — c(l + z) H{z')~^dz' with 
H{z') = Ho[{l + z'fQmM + ^amY^'^. Throughout we adopt 
the cosmological parameters as fim.o ~ 0.27, Qa.o ~ 0.73, 
and Hq = 71 km s~^Mpc~^ (Komatsu et al. 2010). 

Firstly, for the GRB rate R{z), it can in principle be 
connected to the cosmic SFR p*{z), since in the coUapsar 
model the formation of each GRB just indicates the death 
of a short-lived massive star. For relatively low redshifts {z < 
4), the SFR can be expressed approximately by (Hopkins & 
Beacom 2006) 



p*iz) oc 



(1+^) 
(1+^) 



3.44 
0.26 



Z < 0.97, 

0.97 < 2 < 4, 



(5) 



with /6*(0) — 0.02 MQyr~^Mpc~'^, whereas the star forma- 
tion history above 2 ~ 4 is unclear so far. So 12 GRBs with 
a > 4 (the data in region C in Figure 1) are excluded from 
our statistics. Secondly, for the GRB LF, two representative 
forms are usually assumed in literature as a BPL 



(t){L) oc 



L < Lt, 
L > Lt. 



(6) 



and a SPLEC 



^ The average fluxes are calculated from 5'[i5_i50keV] 7^90 where 
'S'[i5— isokcV] is the observed fluences. 



(j,{L) oc 



-Lp/L 



(7) 



Then the expected number of GRBs with redshift 21 < 2 < 
22 and luminosity Li < L < L2 can be calculated by 



max[Li,itiJz)] 



il + ,)^^iL)p,{z)dL^ 



(8) 



where the extra evolving factor (1 + 2)^ is introduced by 
considering that (i) the connection between the GRB rate 
and the SFR could be not in a trivial way and (ii) the LF 
could evolve with redshift. Corresponding to different selec- 
tion criterions and bin methods for different GRB samples, 
the specific form of the above equation (e.g., the sequence 
and the range of the integrals) should be changed, see Equa- 
tions (10), (11), (13), (15), and (18). 

The luminosity threshold invoked in Equation (|8]) can 
be given by 



Lth(2) =47rdi(2)2Pthfc(2), 

wherefc(.) = /,tv^^S(£')i?'rfi5'//, 



150{l-l-z)keV 



(9) 

S{E')E'dE' 



IlkcV "\ J I Jl5{l+z)kcV 

(the primes represent rest-frame energy) converts the ob- 
served flux in the BAT energy band 15 — 150 keV into the 
bolometric flux in the rest-frame 1 — 10'' keV. The observed 
photon number spectrum S(E) can be well expressed by 
the empirical Band function (Band et al. 1993), more 
simply, a broken power law. The value of k varies from 
5.4 to 2.1 as the redshift increases from to 10, by taking 
the rest-frame peak energy as E'^ ^ 200 keV (the most 
frequent value in the Butler et al.'s database) and the 
photon indices prior and post the break energy as 1 and 
2.25 , respectively (Preece et al. 2000). On the other 
hand, unfortunately, a precise description for Pth is nearly 
impossible, since the trigger of the BAT is very complicated 
and, especially for GRBs with redshift, the actual threshold 
is determined by the combination of the BAT and other 
related telescopes. Instead of an abrupt cutoff at Pth, a 
realistic situation could be that the detection efficiency 
starts to remarkably decrease at a certain fiux and ap- 
proaches to be zero with decreasing fluxes. Therefore, in 
the following calculations, a relatively high value for Pth is 
taken as 5 x 10~^ erg s~'^cm~^ and, correspondingly, 12 
GRBs below the selected threshold are further excluded (as 
shown in Figure 1). Strictly speaking, the selected Pth is 
not but higher than the true sensitivity. This can make us 
avoiding the complex arising from the trigger probability. 
The availability of the selected Pth will be checked by a 
fltting to the log A'^ — log P distribution in section 4.3. 



2.3 Observational luminosity distribution 

As the main objective we want to account for, the luminos- 
ity distribution of the selected 101 GRBs (solid asterisks in 
Figure 1) is presented in Figure 3, where an obvious break 
appears at 3 x 10^'erg s~'. Such a break may reflect an 
intrinsic break in LF or just be caused by the selection ef- 
fects arising from the BAT and also other related telescopes, 
which is what we want to clarify in this paper. In order to 
avoid the consideration of the flux-dependence of the red- 
shift measurements in our analyses, we count an effective 
number as i9j^(P) for each GRB with flux P whose redshift 
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Figure 3. Luminosity distribution of the selected 101 GRBs, 
where the horizon error bars correspond to the bin width. Open 
and solid circles correspond to pre- and post-corrected distribu- 
tions, respectively, and the solid line gives a power-law fitting for 
the corrected high-luminosity distribution. The equivalent num- 
ber of the corrected sample is about 319. 

is measured. As a result, an effective GRB sample of a num- 
ber of about 319 is derived from the 101 GRBs. Since darker 
GRBs have higher weight in the statistics, the corrected lu- 
minosity distribution becomes steeper, especially above the 
break luminosity. A good power-law fitting as A'" cx to 
the high-luminosity distribution indicates 1^2 = 1.98 in the 
BPL model and v — 1.98 in the SPLEC model, in view of 
the probable unimportance of most observational selection 
effects at high-luminosity range. In the following Figures 4, 
6, 7, and 8, the observational data are all with the same 
number correction. 



3 EVOLUTION EFFECT 

Before we use the observational luminosity distribution to 
constrain the GRB LP with equation ((8]), we should de- 
termine the evolution parameter A in advance. Following 
Yiiksel et al. (2008) and Kistler et al. (2009), the value of A 
can be constrained by fitting the observational cumulative 
redshift distribution of GRBs with relatively high luminosi- 
ties {L > Lcut). The cut luminosity Lcut is chosen to be 
equal to or higher than the threshold at the highest redshift 
of the sample (here Zmax = 4), so that, in the corresponding 
theoretical calculation, the integral of the LF can be treated 
as a constant coefficient no matter the specific form of the 
LF, i.e., 

iV<?(x r""'%(L)dL r(l + ^rP.(-')^^. (10) 

Due to the limited size of the sample, the observational red- 
shift distribution actually is slightly dependent on the se- 
lected Lcut. So we take Lcut = Lth(4) = 2 x 10^^ erg s~^ to 
reduce the statistical uncertainty as much as possible. Con- 
sequently, 63 GRBs (the data in region A in Figure 1) are 
selected. A comparison between the model and the observa- 
tion is presented in the upper panel of Figure 4, which shows 
that the non-evolution case (A — 0) can be definitely ruled 
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Figure 4. Normalized cumulative number of the 63 GRBs with 
z < 4 and L > 2 X 10^^ erg s~^ as a function of redshift (his- 
togram). Upper panel: Fittings to the redshift distribution with 
Equation l|10p . where the solid and dashed lines correspond to 
A = 1.93 and A = 0, respectively. Lower panel: Fittings to 
the redshift distribution with Equation lllll l. where the solid and 
dashed lines correspond to ^max = 0.72Zq, S = 1.56 (best fit) 
and Zmax = O.IZq, <5 = 0, respectively. 

out, as found before (e.g., Salvaterra & Chincarini 2007; Sal- 
vaterra et al. 2009; Kistler et al. 2008, 2009). The best fitting 
to the observation gives A = i.gfl 

Then an interesting question arises as where such an 
evolution comes from. As found by MacFadyen & Woosley 
(1999) and Yoon et al. (2006), the formation of the black hole 
(or neutron star) during the collapse can drive a GRB event 
only if the coUapsar has high angular momentum. In order to 
avoid strong stellar winds loosing angular momentum, GRB 
progenitors are required to be in low-metallicity evironments 
(Woosley & Heger 2006; Yoon et al. 2006). This theoreti- 
cal metallicity requirement is widely favored by the estima- 
tions of the metallicities of long GRB hosts (e.g. Chen et al. 
2005; Gorosabel et al. 2005; Starling et al. 2005). Therefore, 
it is suggested that the observationally required evolution 
could be mainly due to the cosmic evolution of metallic- 
ity. Specifically, as derived by Langer & Norman (2006), the 
fraction belonging to metallicity below Zmax can be calcu- 
lated by *<z„,.j2) = t[0.84, {Z^,JZq)HO°-^-] /r(0.84), 

* If we do not correct the GRB number by the factor i)J^{P), we 
can get A = 1.44, which is consistent with the results in Kistler 
et al. (2008, 2009). 
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Figure 5. x^-probability distribution of the fittings with Equa- 
tion mil l to the redshift distribution of the 63 GRBs with z < 4 
and L > 2 X lO^^erg s~^. The best-fit parameters 5 = 1.56 and 
^max = 0.72Zq is labeled by the cross. The hatched region repre- 
sents the theoretical metallicity range expected by some coUapsar 
models. 



where ^max is the maximum metaUicity available for GRB 
progenitors and F and F are the upper incomplete and com- 
plete gamma functions. Following this consideration, Equa- 
tion (|10|l becomes 



NT? 



(t>{L)dL 



(1 + /)"* 



<z„ 



dVjz') 
1 + z' 



Varying the parameters Zmax and 5, we refit the observa- 
tional redshift distribution shown in Figure 4 and present 
the distribution of the x^-probabilities of the fittings in Fig- 
ure 5. At first sight, the best-fit parameters .^max = 0.72Zq 
and S — 1.56 may indicate that there is a significant ex- 
tra evolution other than the metallicity evolution. However, 
the long and narrow contours shown in Figure 5 robustly 
demonstrate that the parameters Zmax and S are actually 
strongly degenerate and, moreover, the specific values of the 
best-fit parameters are probably sensitive to the selection of 
the observational sample. Therefore, instead of paying at- 
tention to the best-fit parameters, we treat the relationship 
between the two parameters exhibited by the contours as a 
more valuable result, very roughly, 5 ~ 2.4(^„iax/^0 —0.06). 
Anyway, an independent constraint on these two parameters 
is demanded in order to reduce the parameter degeneracy. 

For example, a theoretical constraint on metallicity can 
be invoked (e.g., Campisi et al. 2010). As proposed by 
Woosley & Heger (2006) and Yoon et al. (2006), the max- 
imum metallicity available for GRB progenitors is likely 
to be within ~ [0.1 — 0.3]^©. As shown in Figure 5, for 
^max < 0.3^0, the value of S would not be higher than 
0.8 with 99.7% confidence. Especially for Zmax < O.IZq, 
the value of 6 approaches to be zero. The fitting to the ob- 
servation with Zmax = O.IZq and 5 = is shown in the 
lower panel of Figure 4 in comparison with the fitting with 
•^max = O.72Z0 and S = 1.56. As can be seen, the difference 
between these two fittings is not very significant. Therefore, 
an extra evolution other than the metallicity evolution may 
be not inevitable if Zmax is indeed very low. 

In the following calculations, we take the best-fit param- 
eters Zmax = 0.72 and 5 — 1.56 just for a good description 



for the evolution effect. Constraints on the LF actually can 
not be significantly affected by the variation of Zmax and S as 
long as they satisfy the required relationship. On the other 
hand, for simplicity, we will ascribe the possible extra evo- 
lution to some unknown factors in the connection between 
the GRB rate and the SFR, i.e.. 



Riz) = Cfl(H-z)**< 



Xz)p*iz), 



(12) 



where the proportional coefficient Cr will be determined 
in Section 4.4. In other words, the LF will be taken to be 
non-evolving in this paper. 



4 LUMINOSITY FUNCTION 
4.1 A preliminary constraint 

With given Zmax and 5, we can constrain the unknown LF 
by fitting the observational luminosity distribution of the 
101 GRBs by 

/ ^(L)R{z)-^dL, (13) 

Li Jo J- + Z 



which gives the expected number in each luminosity bin 
Li < L < L2- The maximum redshift zm{L) as a function 
of luminosity can be solved fron|3 



(11) 



c{l + z) 
Hq 

i ^ 



1 



^(l-f z)^rim,0 + f^A.O 
1/2 



zdz 



\47vPthkJ 



(14) 



With fixed U2 = 1.98 in the BPL model and u = 1.98 in the 
SPLEC model and minimizing the x'^-statistic of the fittings, 
we obtain the best-fit parameters as Lb — 2.5 x lO^^erg s~^ 
and 1/1 = 1.67 for the BPL model and Lp = 2.5 x lO^'^erg s~^ 
for the SPLEC model. As shown in Figure 6, the fitting with 
a BPL LF seems much better than the one with a SPLEC 
LF. This impels us to favor the BPL model. However, the 
apparent oscillation of the observational data, which can not 
be explained by both models, still demands a much more 
elaborate fitting. 

4.2 Redshift-desert effect 

With the preliminary LFs derived above, we can give some 
model-predicted cumulative redshift distributions by 



(t>{L)R{z')dL 



JLa,M 



dV{z') 
1 + z' ' 



(15) 



which are shown in the upper panel of Figure 7 in compari- 
son with the observational one of the 101 GRBs. Obviously, 
the observational numbers at the middle redshifts are much 
less than the ones predicted by both the models. Such a 
remarkable dip in the observational redshift distribution is 
probably, at least partly, related to the so-called 'redshift- 
desert' effect, which is ignored in the above analyses. As 

^ With an approximate expression for luminosity dis- 
tance as di{z) Si W'^^ + z{y/l + z — 1), the max- 
imum redshift can be approximatively calculated by 
ZM ~ 5 {y/l + 4Hodi^M/3c + 2Hodi_M/3c - l) . 
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Figure 6. The best fitting to the observational GRB luminos- 
ity distribution with equation M'il . The dotted and solid lines 
represent intrinsic and observable distributions, respectively. The 
parameters are Li, = 2.5 X lO^^erg , vi = 1.67, and U2 = 1.98 
for the BPL model (upper) and Lp = 2.5 X lO'^^crg s~^ and 
ly = 1.98 for the SPLEC model (lower). 
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Figure 7. Model-predicted redshift distributions with (upper) 
and without (lower) redshift-desert effect in comparison with the 
observational distribution, where the adopted LPs are the same 
to Figure 6. The parameters for the redshift-doscrt effect read 
H = 1.80 and a = 0.79 for the BPL model and fi = 1.63 and 
a = 0.61 for the SPLEC model. 



qualitatively analyzed by Fiore et al. (2007), it could be dif- 
ficult to measure redshifts within the range 1.1 < z < 2.1, 
since at z > 1.1 some strong observable emission or ab- 
sorption lines are shifted outside the typical interval cov- 
ered by optical spectrometers (3800 — 8OOOA) while Lyman- 
Q enters the range at z ~ 2.1. In this paper we do not 
try to give a theoretical description for the redshift-desert 
effect, which must involve many physical and technical is- 
sues. We also notice that a same significant dip actually can 
not be found in the redshift distribution of only GRBs with 
L > 2 X 10^^ erg s~^, as shown in Figure 4. Hence, we sus- 
pect that the redshift-desert effect may mainly influence the 
redshift measurements of relatively dark GRBs, which can 
also be implied by the luminosity distribution. 

Therefore, we show the observational redshift distribu- 
tion of only 38 GRBs with z < 4 and L < 2 x lO^^erg s'^ 
(the data in region B in Figure 1) in the lower panel of Fig- 
ure 7. Meanwhile, for relatively dark GRBs, we tentatively 
suggest a Gaussian function 



77,(2) = 1 — exp 



{z-^r 



(16) 



to phenomenologically describe the redshift-dependence of 
the redshift measurements. In contrast, for sufficiently bright 



GRBs, we take 
'7.(2) = 1- 



(17) 



Fittings to the distribution of the 38 GRBs give the best-fit 
parameters as /i = 1.80 and u = 0.79 for the BPL model 
and /i = 1.63 and a = 0.61 for the SPLEC model, which 
are basically consistent with the theoretical expectation of 
the redshift-desert effect. With these phenomenological ex- 
pressions of riz{z), we refit the redshift distribution of the 
101 GRBs, which is also shown in the lower panel of Fig- 
ure 7. As can be seen, the fittings are greatly improved, as 
the observational dip in the redshift distribution is produced 
naturally, especially in the BPL model. 



4.3 Final results 

Combining Equations fTS}, (|16p . and (|17p . we refit the lu- 
minosity distribution of the 101 GRBs and find that Li, = 
2.5 X lO^^erg s"^ and ui = 1.72 for the BPL model and 
Lp = 2.5 X 10"^ erg s"^ for the SPLEC model. Strictly, we 
should use these new LFs to re-constrain the redshift-desert 
effect and go recycling until reaching a certain precision. For 
simplicity, we stop here because the obtained new values of 
the parameters are only slightly different from the prelimi- 
nary ones. With these new parameters. Figure 8 shows that 
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Figure 8. The same to Figure 6 but with a redshift-desert ef- 
fect derived from Figure 7. The LF parameters are Lj = 2.5 X 
lO^^erg s-\ i/i = 1.72, and V2 = 1-98 for the BPL model and 
Lp = 2.5 X lO^'^erg s"! and u = 1.98 for the SPLEC modeh 



the BPL model complies with the observation successfully, 
whereas the SPLEC model still predicts some remarkable 
excesses around ~ 10^° - lO^^erg s~\ Therefore, we prefer 
to conclude that the GRB LF could be a BPL. 

Finally, with the derived BPL LF, we give a model- 
predicted cumulative flux distribution by 



(j){L)R[z)-&^[P')r)^[z)dL 



dV{z) 

1 + 2 



(18) 



where Lp = AirdikP, P' = L/Andfk, and the redshift de- 
tection probability '&z{P')riz{z) has been determined above. 
As shown in Figure 9, the consistency between the theoret- 
ical and observational flux distributions indicates that our 
choice of the luminosity threshold is basically reasonable, 
i.e., the trigger probability above Pth by the BAT can be 
afhrmed to nearly constant. 




-8.0 -7.5 -7.0 -6.5 
log (P / erg s ^cm ^) 



5.5 



Figure 9. Comparison of the modcl-prcdictcd flux distribution in 
the BPL LF model with the observational one (without number 
correction) of the 101 GRBs. The model parameters are the same 
to Figure 8. 
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.dV{z) 



1 + z 
(19) 

where (An/47r) ~ 0.1 is the field of view of the BAT, T ~ 5 
yr is the observational period, fb ~ 0.01 is the beaming de- 

gree of the GRB outflow, and ~ {vi — V) { —f^ j 77^ 
is the normalization coefficient of the BPL LF with Lmin 
being an assumed minimum luminosity for the GRBs. 

For & = 1.56, Zniax = 0.72^0, = 1.72, U2 = 1.98, 
Lb = 2.5 X lO^^erg s~\ and N^53,25,53.75] = 1-9, the propor- 
tional coefficient in the GRB rate can be constrained to 



Cb = 4 X 10' 



W-'^Lb 



h 

0.01 



M, 



' 



(20) 



Gpc "^yr ^ and an ob- 



which yields a overall local GRB rate as 
Ci?*<z„,jO)p.(0) = 36(/6/0.01)-^ 
served local GRB rate as fbR{0) = 0.36 Gpc"^yr"\ This 
rate is basically consistent with the previous results (e.g., 
Schmidt 2001; Guetta et al. 2004, 2005; Liang et al. 2007; 
Wanderman & Piran 2010). The value of Cr also implies 
that, besides the metallicity requirement, GRB progeni- 
tors may also have some other particularities. For exam- 
ple, as widely accepted, only massive Wolf-Rayet stars (e.g., 
> 20Mq) are possible GRB progenitors (MacFadyen & 
Woosley 1999; Larsson et al. 2007). So a small fraction 

arises as /wa = SIom '^^ '^^'^)'^'^ / /o^imq® "^^("^)'^'^ ~ 
2 X 10""^ -^0^1 where •^{m) is the Salpeter initial stellar 
mass function. Additionally, there is still an extra factor of 
~ 10^^ — 10^^ unexplained, which could be related to the 



4.4 GRB rate 

In the above fittings to the luminosity distribution, we nor- 
malize the model-predicted GRB number by the following 



° The normalization is usually estimated with an entire dataset or 
a good statistical point. Here, we select the data at the highest lu- 
minosity L = 10^^'^erg for normalization due to two reasons: 
(i) the data above L = lO^^erg s~^ can be well fitted by a power 
law, which indicates that all the data above L = 10^'^erg s~^ are 
probably good statistical, and (ii) more higher-luminosity GRBs 
may have loss selection effects. 
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particularity of GRB progenitors in their rotation, magnetic 
fields, etc. 



5 CONCLUSION AND DISCUSSION 

The accumulation of Swift observed GRBs gradually makes 
it possible to directly derive a GRB LF from observational 
luminosity distribution, where however two complexities 
must be involved as (i) the evolving connection between 
the GRB rate and the cosmic SFR and (ii) observational 
selection effects. With a phenomenological investigation on 
these two complexities, we constrain and discriminate two 
popular competitive LF models and find that the BPL LF 
model is more favored by the observation. However, in view 
of the approximative description of the selection effects, the 
SPLEC still can not be ruled out absolutely. 

Although the derived values of the parameters ji and 
a are basically consistent with the theoretical expectation 
of the redshift-desert effect, the flux-dependence of the 
redshift-desert effect is still very ambiguous (an abrupt lu- 
minosity boundary as 2 x fO^^erg s^^ is adopted in our anal- 
yses). More seriously, the flux- and redshift-dependences of 
the redshift measurements actually must be coupled with 
each other, but in our analyses the functions '&z{P) and 7)z(z) 
are considered independently for simplicity. This may lead 
to an overestimation of the redshift selection effect. There- 
fore, some more elaborate theoretical considerations on the 
redshift measurements are demanded. On the other hand, of 
course, a more detailed analysis on the observational results 
would be helpful, specially towards to every kind of redshift 
measurement methods. 

Finally, our investigation on the evolution effect shows 
that, if the matallicity of GRB progenitors is lower than 
~ O.IZq as expected by some coUapsar models, there may be 
no extra evolution effect (i.e., i5 ~ 0) other than the metal- 
licity evolution. Alternatively, if we remove the theoretical 
metallicity requirement, a stronger extra evolution would be 
required for higher metallicties. In the latter case, the extra 
evolution could indicate some other evolutions in the GRB 
rate or indicate an evolving LF which is not considered in 
this paper. To discriminate these two possibilities is difficult 
but interesting. It will be helpful to separate the GRB sam- 
ple into few subsamples with different redshift ranges and fit 
the luminosity distribution of each subsamples one by one 
with an evolving LF. Such a further work can be made as 
the GRB sample becomes large enough. 
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